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COMPUTER  PROGRAM  FOR  THE  REDUCTION 
AND  PRELIMINARY  ANALYSES 
OF  RUNOFF  DATA\_f 

^V  i     p-  °  ° 

by  \  D.  A.1woolhiser2  and  K.  E.lsaxtonl 

INTRODUCTION 

Electronic  digital  computers  have  been  used  for  the  reduction  and  analyses  of  runoff  data 
for  several  years.  Their  speed  and  accuracy  have  made  them  extremely  valuable  for  this  purpose. 

Most  of  the  available  programs  have  been  written  for  computers  having  smaller  storage 
capacities  than  those  currently  available.  Because  of  their  limited  storage  capacity,  two  runs 
are  usually  required  with  the  output  of  the  basic  data-reduction  program  serving  as  input  for  the 
analyses  program.  In  many  cases,  the  input-output  operations  require  more  time  than  the 
computations. 

This  program  performs  the  reduction  and  preliminary  analyses  of  runoff  data  with  a  single 
pass  of  the  data.  From  the  time-stage  data  of  the  flow  events  and  the  station's  rating  infor- 
mation, the  program  computes  flow  rates,  flow  volumes,  flow  durations  and  volumes  above 
selected  rates,  and  maximum  amounts  of  flow  for  selected  time  intervals. 

The  program  was  written  to  compute  runoff  from  agricultural  watersheds  as  measured 
by  weirs  and  pond  outlet  structures.  However,  it  is  sufficiently  general  that  it  may  be  applied 
to  nearly  all  types  of  flow-measuring  devices  having  a  known  stage-discharge  relation.  Many  con- 
stants within  the  analyses  portion  of  the  program  may  be  easily  modified  to  allow  for  specific 
applications. 

The  program  was  written  in  FORTRAN  for  use  on  a  Control  Data  Corporation  1604  com- 
puter^; however,  it  is  adaptable  to  most  other  computers.  The  authors  have  successfully 
adapted  it,  without  major  programming  changes,  to  an  IBM  1620  having  60,000  units  of  storage. 
This  version  of  the  program  may  be  obtained  from  the  authors. 

Since  the  program  is  generally  adaptable,  its  description  will  be  useful  to  those  engaged 
in  hydrologic  research  involving  field  measurement  of  runoff.  Additional  clarification  of  the 
program  can  be  obtained  from  either  of  the  authors. 


Soil  and  Water  Conservation  Research  Division,  Agricultural  Research  Service,  in  cooperation  with  the  Wis- 
consin Agricultural  Experiment  Station. 

2  Formerly  with  Soil  and  Water  Conservation  Research  Division  at  Columbia,  Mo.;  now    at    Cornell     University, 
Ithaca,  N.Y. 

3  Soil    and    Water   Conservation    Research  Division  at  Columbia,    Mo. 

Trade  names  are  used  in  this  publication  solely  to  provide  specific  information.  Mention  of  a  trade  name 
does  not  constitute  a  guarantee  or  warranty  and  does  not  signify  that  the  product  is  approved  to  the  exclusion  of 
other  comparable  products. 


PROGRAM  FEATURES 


This  program  computes  instantaneous  rates  of  discharge  from  simultaneous  stage  and  time 
readings  and  numerically  integrates  these  rates  with  respect  to  time  to  obtain  runoff  volumes. 
These  volumes  are  summed  by  storms,  days,  months,  and  years.  Flow  volumes  and  durations 
at  discharges  equal  to  or  greater  than  twenty  specified  rates  are  computed,  and  annual  maximum 
volumes  are  determined  for  nine  specified  time  intervals.  All  calculations  and  summaries  for 
reduction  and  analyses  are  accomplished  with  a  single  pass  of  the  input  data. 

Four  optional  methods  of  computing  discharge  rates  are  provided: 

1.  As  a  function  of  stage  alone. 

2.  As  a  function  of  stage  and  ponding. 

3.  As    a   function    of   stage    and  ponding  with  the  rate  of  ponding  corrected  for  rain  falling 
on  the  pond. 

4.  As    a   function    of    stage,    ponding,    rainfall-on-the-pond,    and    a  variable  drainage  area 
caused  by  changes  in  the  surface  area  of  the  pond  as  the  stage  changes. 

A  detailed  description  of  the  computer  program  is  presented  in  the  appendixes.  A  complete 
flow  chart  is  shown  in  appendix  1.  The  flow  chart  has  been  separated  into  14  divisions  according 
to  various  functions.  Some  overlap  is  inevitable,  particularly  in  indexing.  The  primary  function 
of  each  division  is  given  in  appendix  2. 

To  complete  the  program's  description,  the  FORTRAN  statements  are  given  in  appendix  3 
and  the  symbol  definitions  are  given  in  appendix  4.  The  program  statements  are  listed  exactly 
as  they  are  used  in  the  program.  Punching  and  stacking  a  deck  containing  these  statements 
will  produce  a  functional  program. 


INPUT 


The  input  to  the  program  is  read  from  standard  data  processing  cards.  It  consists  princi- 
pally of  the  following: 

1.  Several  cards  of  basic  time  units  and  constants  used  in  computing. 

2.  A  "station  card"  containing  codes  and  information  for  the  particular  type  of  calculations 
to  be  used  for  the  station's  data. 

3.  Tabular  values  defining  the  station's  stage-discharge  relation,  stage-pondage  relation, 
and  stage-area  relation.  The  latter  two  may  be  used  optionally. 

4.  Data  of  either  time-stage  or  time-stage-precipitation  readings  from  recorder  charts. 

5.  Code  cards  for  designating  specific  functions. 

An  example  of  input  to  the  program  is  shown  to  aid  in  understanding  the  program.  Figure  1 
shows  the  instrument  chart  of  time-stage  relation,  and  figure  2  is  a  list  of  the  complete  input 
to  the  program.  This  represents  the  input  for  a  typical  isolated  runoff  event  over  a  weir  on 
an  intermittent  stream  having  a  small  pond  created  by  flow.  Pondage  is  considered;  however, 
rairtfall-on-the-pond  and  a  variable  drainage  area  are  not  considered  since  the  pond  is  small 
in  relation  to  the  drainage  area.  Two  lines  of  input  data,  the  first  and  fifth,  have  been  left  unfilled 
to  illustrate  that  the  input  cards  do  not  have  to  be  full;  two  time  increments,  those  underscored, 
have  intentionally  been  made  negative  to  illustrate  the  elimination  of  the  values  in  error. 

To  adequately  present  each  portion  of  the  input,  a  detailed  description  follows. 
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Figure  2. — Computer  input  of  the  example  event. 


Basic  Constants 


The  basic  constants  used  in  the  computations  are  listed  on  five  cards  as  shown  in  table  1. 

Cards  Nos.  1  and  2  contain  the  accumulated  days  of  the  year  at  the  beginning  of  each  month. 
LDAY(M)  is  for  a  non-leap  year  and  LEAP(M)  is  for  a  leap  year. 

Cards  Nos.  3  and  4  contain  the  rates,  CLASS(J),  in  inches  per  hour,  used  as  limiting  rates 
in  the  flow-duration  computations.  Volumes  and  durations  of  flow  occurring  at  rates  equal  to  or 
greater  than  CLASS(J)  are  computed.  These  rates  may  be  varied  as  desired  by  changing  the 
CLASS(J)  values. 

Card  No.  5  contains  the  time  periods,  P(J),  in  hours,  used  when  searching  for  maximum 
amounts  of  runoff  for  specified  time  periods.  These  values  may  also  be  varied  as  desired. 


TABLE  1. —  Basic  constants  used  in  computations 


Accumulated  days  of  the  year  (format  4)1 

LDAY(M)         0       31        59       90       120       151        181        212       243       273       304       334 
LEAP(M)         0       31        60       91        121        152       182       213       244       274       305       335 

Flow-duration  class  rates,  inches  per  hour  (format  151) 

CLASS  (J)2        .0010        .0016        .0025        .0040        .0063        .010        .016        .025        .040        .063 
.100  .16  .25  .40  .63  1.00        1.60        2.50        6.30      10.00 

Time  lengths  for  maximum  amount  search,  hours  (format  3) 
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1  For  format  statements,  see  FORTRAN  program  listing,  appendix  3. 

2  CLASS  (J)  intervals  may  be  modified  by  changing  the  values  in  this  table.  Those  shown  plot 
evenly  on  a  logarithmic  scale. 

3  The  time  lengths  may  be  set  at  any  selected  value. 

Station  Information 

The  following  is  the  format  and  definition  of  the  constants  provided  on  the  station  information 
card.  The  first  line  indicates  the  spacing  across  a  data  card;  the  second  line  represents  the  signs, 
numbers,  and  decimals  as  they  would  be  used;  and  the  third  lists  the  symbol  for  each  station 
constant  (See  appendix  4  for  symbol  definitions.) 

Rating  Table 

Values  of  the  station's  discharge,  in  units  of  cubic  feet  per  second,  are  read  into  and  stored 
in  the  computer.  These  values  are  to  be  taken  at  uniform  intervals  of  stage  from  the  stage  dis- 
charge relation,  with  the  first  value  being  at  the  station's  gage  datum  elevation,  GDE.  They  are 
punched  in  eight-column  fields  and  ten  fields  per  card  (10F8.3).  The  decimal  may  be  punched,  or 
the  machine  will  establish  it  as  three  places  from  the  right  side  of  the  field  if  unpunched.  All 
cards  must  be  filled  except  the  last. 

The  program  applies  a  linear  interpolation  between  tabulated  values  of  discharge  to  compute 
the  discharge  rate  for  stages  between  two  tabulated  values.  Therefore,  the  stage  interval,  QI, 
between  the  discharge  rates  stored  should  be  selected  so  a  linear  interpolation  will  give  the 
desired  accuracy. 


The  amount  of  storage  set  aside  in  the  computer  for  the  rating  table  values,  Q,  is  specified  in 
the  DIMENSION  statement.  This  amount  should  be  large  enough  to  allow  for  the  maximum  number 
of  entries  expected  but  not  excessive.  The  same  is  true  for  storing  the  pondage  table  values,  PC, 
and  the  watershed  area  table  values,  AREA. 


Pondage  Correction  Table 

If  a  pondage  correction  is  to  be  applied,  the  values  are  punched  in  the  same  format  and 
under  the  same  rules  as  the  rating  table.  The  unit  to  be  used  is  cubic  feet  per  second  for  a  rate 
of  change  in  stage  of  one  foot  per  minute.  The  first  value  must  correspond  to  the  pond  datum 
elevation,  PDE,  specified  on  the  station  information  card. 


Watershed  Area  Table 

If  the  permeable  area  of  a  watershed  varies  significantly  with  changes  in  pond  elevation,  a 
variable  drainage  area  table  may  be  necessary.  Values  in  acres  of  watershed  are  punched  in  the 
same  format  and  under  the  same  rules  as  the  rating  table.  The  first  value  must  correspond  to  the 
pond  datum  elevation,  PDE. 


Data  Cards 


The  time- stage  data  are  punched  in  one  of  two  formats  depending  on  whether  rainfall  values 
are  required. 

The  following  shows  the  information  required  and  the  data  card  format  if  no  rainfall  values 
are  needed;  i.e.,  rainfall-on-the-pond  correction  is  not  applied: 
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All  symbols  are  defined  in  appendix  4.  Note  that  four  time- stage  readings  may  be  placed  on 
one  card.  The  card  does  not  need  to  be  full;  however,  all  values  must  be  in  the  farthest  left  field 
available.  A  new  card  is  required  if  ID,  M,  KDAY,  or  KYR  change. 

The  following  is  the  information  required  and  format  to  be  used  if  rainfall  values  are 
needed;  i.e.,  the  rate  of  change  in  stage  for  computing  the  pondage  rate  is  to  be  corrected  for 
rainfall-  on-  the-  pond . 
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Note  that  only  three  time-stage-rainfall  readings  per  card  are  possible  with  this  format. 

It  may  be  more  convenient  and  efficient  to  change  the  input  format  from  that  shown  if  the 
user  prefers  a  different  time  or  stage  representation.  This  change  would  only  require  modifica- 
tions of  the  input  format  and  of  those  statements  that  convert  time  into  decimal  hours  of  the 
year  and  stage  into  head. 

Codes 

The  following  code  cards  and  tabulating  details  are  used  in  connection  with  data  input. 

1.  To  signify  the  end  of  a  station's  data: 

Use  a  data  card  containing  only  the  station's  identifying  number,  ID,  and  the  code  99  in 
the  spaces  provided  for  M.  This  card  causes  the  maximum  volume  search  to  be  made  on 
the  year's  data  just  computed  and  transfers  control  to  the  beginning  of  the  program  for  a 
new  station  card  to  be  read. 

2.  To  add  or  subtract  a  constant  from  the  stage  readings: 

Insert  a  data  card  containing  only  the  station's  ID,  the  code  55  in  the  spaces  provided 
for  M,  and  the  correction  to  be  made,  SHIFT,  in  the  spaces  provided  for  H(l).  This  cor- 
rection value  will  be  added  algebraically  to  all  subsequent  gage  height  readings  until 
another  correction  or  a  zero  is  inserted  by  a  similar  card.  This  value  is  set  to  zero  at 
the  beginning  of  the  calculations  for  each  station. 

3.  To  set  the  storm  total  and  rates  to  zero: 

Place  a  negative  sign  on  an  hour  reading,  TH(I).  This  sets  the  storm  volume  to  zero  and 
the  corrected  discharge  at  this  time  and  the  previous  time  to  zero.  A  large  accumulation 
of  flow  may  be  calculated  for  a  time  of  no  flow  in  computations  involving  reservoirs 
because  the  pondage  rate  is  not  exactly  zero.  This  occurs  because  the  pondage  rate  is 
calculated  from  three  consecutive  stage  values  to  arrive  at  a  rate  for  the  central  value. 
The  pondage  rate  will  not  be  zero  unless  all  of  these  values  are  equal  or  the  rates  of 
change  in  stage  before  and  after  the  central  value  are  equal  and  with  opposite  signs.  Of 
course,  this  code  should  not  be  used  unless  the  discharge  rate,  inflow  to  a  reservoir, 
should  in  fact  be  zero. 

4.  To  stop  the  program: 
Insert  a  blank  card. 

Special  Notes 

Due  to  the  method  of  calculation  used  in  this  program,  the  last  time-stage  reading  of  a  sta- 
tion's data  is  not  printed  when  changing  stations  or  at  the  end  of  the  data;  i.e.,  when  the  code 
M=99  is  used.  To  avoid  losing  this  last  reading,  a  dummy  reading  is  added  after  the  last  actual 
reading.  It  should  have  a  time  a  few  minutes  later  and  the  same  gage  height  as  the  last  actual 
reading. 

The  first  time  increment  of  a  new  station  year  is  artificial  and  has  no  significance. 

No  provision  has  been  made  for  estimated  readings,  but  a  code  such  as  a  negative  minute 
reading  could  easily  be  written  into  the  routine  if  desired.  A  provision  to  handle  estimated 
volumes  could  also  be  added. 

If  a  time  difference  between  two  successive  readings  becomes  negative  because  of  a  tabulat- 
ing or  punching  error,  the  program  will  eliminate  readings  until  the  time  increment  again 
becomes  positive.  When  this  occurs,  an  error  statement  is  printed  and  the  computations  are 
continued. 


Order  of  Input  Cards 

The  input  cards  must  be  arranged  as  follows: 

1.  Program  cards  (FORTRAN  statements). 

2.  Cards  that  may  be  specified  by  the  computer  residence. 

3.  Basic  data  cards: 

LDAY(M)  (one  card) 

LEAP  (M)  (one  card) 

CLASS  (J)  (two  cards) 

P  (J)  (one  card) 

4.  Station  card. 

5.  Discharge  rating  table. 

6.  Pondage  rating  table  (if  required). 

7.  Variable  drainage  area  table  (if  required). 

8.  Data    (including  code  card  M=55  if  required).  This  may  include  data  for  more  than  one 
year  for  the  same  station. 

9.  Code  card  M  =  99. 

10.    The  next  station  card  or  a  blank  card  to  stop  the  program. 
Repeat  card  groups  4  through  9  for  each  station's  data  to  be  computed. 


METHODS  OF  CALCULATION 

The  program  performs  three  principal  calculations.  They  are:  (1)  computing  flow  rates  and 
integrating  them  with  respect  to  time  to  obtain  runoff  volumes;  (2)  determining  duration  and 
volume  of  flow  above  specified  rates  of  flow,  CLASS  (J);  and  (3)  searching  for  maximum  volumes 
of  flow  occurring  during  specified  lengths  of  time,  P  (J).  The  following  is  a  brief  discussion  of 
each  of  the  calculations. 

Flow  Rates  and  Volumes 

To  determine  flow  volumes,  V,  the  program  first  determines  instantaneous  flow  rates  for 
each  stage  reading  by  referring  to  the  station's  rating  table.  Incremental  volumes  are  obtained 
by  averaging  two  successive  rates  and  multiplying  by  the  increment  of  time  between  them,5  thus: 

Av  =  qo  +  qob  xAti, 

as  shown  in  figure  3. 

If  ponding  is  considered,  the  rate  of  flow  is  adjusted  to  account  for  the  amount  ponded  by 
computing  the  rate  of  change  in  stage  at  HC1  by  the  relation: 


RATE  =  |"H£1^HC0       HC2  -  HCll 

L       At,  At2        J       ' 


Atx  At2 

This  approximates  the  slope  of  the  tangent  on  the  time-stage  curve  at  HC1  by  averaging  the 
point-to-point  slopes  on  either  side  of  HCl.If  rainfall-on-the-pond  is  considered,  the  stage  read- 
ings are  reduced  by  the  amount  of  rainfall  before  RATE  is  calculated. 


5  A  more  detailed  description  of  this  method  of  computing   flow    volumes    may  be  found  on    page  153   of  USDA 
Agriculture  Handbook  No.  224,  "Field  Manual  for  Research  in  Agricultural  Hydrology." 
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QCORR 
or   QOB 


l-At,  — u-aIH 


TIME,  t 

Figure  3. — Calculation  of  flow  volumes. 

This  rate  of  change  in  stage,  RATE,  is  then  multiplied  by  the  value  in  the  pondage  table 
corresponding  to  stage  HC1  to  obtain  a  pondage  correction  flow  rate,  QP.  This  value  is  added 
algebraically  to  the  observed  rate,  QOB,  to  obtain  a  corrected  rate,  QCORR.  Therefore: 


and 


QCORR  =  QOB  +  QP 


Av  =  QCORR  +QOxA(i 


These  increments  of  volume  are  then  accumulated  to  provide  the  storm,  daily,  monthly,  and 
annual  totals.  The  accumulated  annual  volumes  and  their  corresponding  times,  in  hours  from  the 
beginning  of  the  year,  are  stored  for  the  entire  year  to  be  used  in  the  search  for  maximum 
volumes  for  specified  time  intervals. 


Flow  Durations  and  Volumes 

After  each  incremental  volume  is  computed,  the  portion  of  the  volume  AV  and  its  correspond- 
ing time  A  t^  that  occur  above  a  flow  rate  CLASS  (J)  are  accumulated  as  VOLD  (J)  and  DUR  (J). 
Therefore,  by  using  several  flow  rates,  CLASS  (J),  as  specified  in  the  basic  data  cards,  a  flow- 
duration  series  is  computed. 

Consider  the  hydrograph  shown  in  figure  4.  The  volume  of  flow  at  rates  equal  to  or  greater 
than  the  rate  CLASS  (J)  is  shaded.  A  straight- line  variation  of  discharge  between  tabulated  points 
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TIME,     t 

Figure  4.—  Hydrograph  for  determining  flow-duration  values. 

is  assumed  and  the  incremental  volumes  and  durations  are  computed  on  this  basis.  For  example, 
the  equation  for  the  incremental  duration  in  triangle  (1)  is: 


*  ,     ^        [Q3  -  CLASS  J  1  x  (t3  -  to) 
Length  (ac)  =   ^ i £1 

Q3-  Q2 


and  the  incremental  volume  is: 

Area.l)    =  fo  -  CLASS  (J)]2  ,  (I3  -  t2) 

2(Q3  "  Q2> 

Maximum  Volumes 

Annual  maximum  volumes  of  runoff  occurring  within  the  selected  time  intervals  P  (J)  are 
determined  by  a  search  of  the  previously  calculated  data.  Two  indexed  variables,  V(L)  and  T(L), 
are  utilized.  V(L)  is  the  accumulated  runoff  in  inches,  and  T(L)  is  its  corresponding  time  in 
hours  from  the  beginning  of  the  calendar  year. 

The  limit  of  storage  in  the  computer  available  for  storing  these  values  is  set  in  the  DIMEN- 
SION statement.  Since  one  unit  of  storage  for  each  variable  is  required  for  each  input  reading, 
an  estimate  can  be  made  of  the  maximum  number  needed.  This  represents  a  sizable  amount  of 
storage  on  the  smaller  computers;  therefore,  it  may  be  necessary  to  store  fewer  values,  search 
more  often,  and  combine  the  results. 

A  mesh  of  time  length  P  (J)  is  set  up  to  make  this  search.  The  mesh  is  indeed  at  J  =  1,  which 
sets  its  length  equal  to  the  first  value  specified  in  the  basic  constants.  This  mesh  is  then  shifted 
through  the  stored  data  to  obtain  the  maximum  volume  accumulated  during  this  length  of  time 
and  the  time  of  the  beginning  of  this  maximum  amount.  The  mesh  is  reindexed  and  the  procedure 
repeated  for  each  successive  length  of  time. 
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For  making  this  search,  LL  indexes  the  leading  edge  of  this  mesh  and  LT  the  trailing  edge. 
These  indices  are  assigned  values  equal  to  the  index  of  a  time  value.  The  time  value  used  is  the 
one  at  which  the  mesh  edge  occurs  if  the  mesh  edge  is  coincidental  with  a  time  value,  or  the  one 
immediately  following  if  the  mesh  edge  occurs  between  two  time  values.  This  substitution  was 
made  in  the  following  three  equations  when  incorporating  them  into  the  program. 

There  are  three  possible  positions  of  the  search  mesh  with  respect  to  the  stored  values  of 
time  as  it  is  shifted  through  the  data: 

Case  A. — The  leading  and  trailing  edges  exactly  coincide  with  the  stored  time  values  T^  as 
shown  in  figure  5.  Thus  no  interpolation  is  required  and  the  volume  Vj  in  the  time  period  P  (J) 
is  given  by: 

Vj  =  vL+N-vL 


Case  B. —  The  leading  edge  falls  between  two  time  values  as  shown  in  figure  6.  An  interpola- 
tion is  required  at  the  leading  edge  of  the  mesh  and  the  volume  Vj  is  given  by: 


V     =  V  + 

J         L+N-l  + 


P(J)  -  T 

L+N-l 


T  -  T 

L+N         L+N 


-1       \ 


V         -  V  \-v 

L+N         L+N-l)       L 


Case  C. —  The    trailing   edge   falls  between  two  time  values,  as  shown  in  figure  7.  Here  an 
interpolation  is  required  at  the  trailing  edge  of  the  mesh;  thus,  the  volume  Vj  is  given  by: 


T  -  T     -  P(J) 

V     =  V  -  V      + x  /  V  -  V 

J         L+N         L  TL+N"TL  '       +1 


(vl+i"vl) 


> 

(TL  +  N'VL+N 

y^ 

o 

_l 

t 

(TL+2,VL+2)___ 

Li. 

Q 

LU 

_l 
=> 

V( 

1 

J) 

x^ui.W 

O 
<) 

D(    1  ^ 

< 

^ 

r\J  ) 

/    LT=L 

LL=L+N 

TIME,  T 


(T        V        ) 

1  'l+n,vl+n' 


o 


Q 

ixl 

I- 
< 


o 
o 

< 


LT=L 


LL= L+N 


TIME,  T 

Figure  6.— Maximum  volume  search  with  leading  edge  interpolation. 


o 

_l 
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T 

V(J) 


(TL+|,VU|) 


3 

o 
o 
< 


(TL,VL) 


^TL+N'VL+N  ^ 


(TL+(N-I)»VL+(N-Ip 


LL  =  L  +  N 


TIME  ,   T 


Figure  7. — Maximum  volume  search  with  trailing  edge  interpolation. 
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OUTPUT 


The  output  of  the  program  includes  a  line  of  data  for  each  time-stage  reading.  Figure  8  shows 
an  example  of  output  obtained  from  the  input  example  shown  in  figure  2. 
Included  in  each  line  are: 

1.  An  identifying  number. 

2.  The  month,  day,  year,  and  time  of  the  reading. 

3.  The  time  increment  between  readings — hours. 

4.  The  corrected  stage — feet. 

5.  The  rate  of  change  in  stage — feet  per  minute. 

6.  The  pondage  correction — cubic  feet  per  second. 

7.  The  observed  discharge — cubic  feet  per  second. 

8.  The  observed  discharge  corrected  for  pondage — cubic  feet  per  second  and  inches  per  hour. 

9.  Accumulated  storm,  daily,  monthly,  and  annual  flow  volumes — inches. 

The  error  statements  indicate  where  the  difference  between  two  time  readings  became  negative 
and  the  latter  reading  was  omitted. 

At  the  end  of  the  year's  data  or  when  the  code  card  M=99  is  inserted  to  signify  the  end  of  a 
station's  data,  the  maximum  volumes  occurring  within  the  selected  time  intervals  P  (J)  are 
printed  as  inches  of  runoff.  The  time  intervals  used  in  the  example  shown  were  arbitrarily 
selected. 

Flow-duration  information  is  printed  after  the  maximum  volumes.  This  includes  the  times, 
in  hours,  and  the  volumes,  in  inches,  that  the  flow  was  equal  to  or  greater  than  the  twenty 
selected  discharge  rates. 
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ID  MO  DAY  YR  HOUR  KEN  DELTTIME  HEAD 
.0167 
.0500 
.0333 
.0333 
.0167 
.0333 
.0333 
.0333 
.0333 
.0667 
.0500 
.0500 
.0833 
NEGATIVE  DELT  TIME  INCURRED  Ai 
152  7  2  60  10.  9.  .0667 
152  7  2  60  10.  22.  .2167 
NEGATIVE  DELT  TIME  INCURRED  AI 
152  7  2  60  10.  30.  .1333 
152  7  2  60  11.  0.  .5000 
152  7  2  60  11.  14.  .2333 
152   7   2  60  11.  16.     .0333 


152 

7 

2 

60 

9. 

34. 

152 

7 

2 

60 

9. 

37. 

152 

7 

2 

60 

9. 

39. 

152 

7 

2 

60 

9. 

41. 

152 

7 

2 

60 

9. 

42. 

152 

7 

2 

60 

9. 

44. 

152 

7 

2 

60 

9. 

46. 

152 

7 

2 

60 

9. 

48. 

152 

7 

2 

60 

9. 

50. 

152 

7 

2 

60 

9. 

54. 

152 

7 

2 

60 

9. 

57. 

152 

7 

2 

60 

10. 

0. 

152 

7 

2 

60 

10. 

5. 

HEAD 

DHDT  PONDCOR  QOB(CFS) 

CFS 

Q(1PH) 

STORM 

DAILY 

MONTHLY 

ANNUAL 

1.32 

.0000 

.00 

.00 

.00 

.0000 

.0000 

.0000 

.0000 

.0000 

1.36 

.0167 

.23 

.00 

.23 

.0099 

.0002 

.0002 

.0002 

.0002 

1.40 

.0250 

.36 

.00 

.36 

.0158 

.0007 

.0007 

.0007 

.0007 

1.46 

.0350 

.55 

.00 

.55 

.0241 

.0013 

.0013 

.0013 

.0013 

1.50 

.0550 

.92 

.00 

.92 

.0400 

.0019 

.0019 

.0019 

.0019 

1.64 

.0600    1 

.20 

.06 

1 

.26 

.0547 

.0034 

.0034 

.0034 

.0034 

1.74 

.0450    1 

.04 

.22 

1 

.26 

.0546 

.0053 

.0053 

.0053 

.0053 

1.82 

.0375 

.99 

.44 

1 

.43 

.0624 

.0072 

.0072 

.0072 

.0072 

1.89 

.0262 

.80 

.72 

1 

.52 

.0662 

.0094 

.0094 

.0094 

.0094 

1.96 

.0104 

.36 

1.10 

1 

.46 

.0637 

.0137 

.0137 

.0137 

.0137 

1.97 

-.0000 

.00 

1.16 

1 

.16 

.0505 

.0165 

.0165 

.0165 

.0165 

1.96 

-.0057 

.20 

1.10 

.90 

.0392 

.0188 

.0188 

.0188 

.0188 

1.92 

-.0103 

.33 

.87 

.53 

.0233 

.0214 

.0214 

.0214 

.0214 

'ER  NEXT  READING 

1.87 

-.0120 

.35 

.63 

.28 

.0120 

.0226 

.0226 

.0226 

.0226 

1.72 

-.0095 

.21 

.17 

.00 

.0000 

.0239 

.0239 

.0239 

.0239 

■ER  NEXT  READING 

1.66 

-.0059 

.12 

.08 

.00 

.0000 

.0239 

.0239 

.0239 

.0239 

1.53 

-.0032 

.06 

.00 

.00 

.0000 

.0239 

.0239 

.0239 

.0239 

1.50 

-.0011 

.02 

.00 

.00 

.0000 

.0239 

.0239 

.0239 

.0239 

1.50 

.0000 

.00 

.00 

.00 

.0000 

.0239 

.0239 

.0239 

.0239 

ANNUAL  FLOW  SUMMARY   ID=  152   YEAR=  60 


MAXIMUM  FLOWS 


DURATION  OF  FLOW  RATES 


ID  INTER 

MO 

DAY 

HOUR 

VOLUME 

152   1. 

7 

2 

9.57 

.024 

152   2. 

0 

0 

.00 

.000 

152   6. 

0 

0 

.00 

.000 

152   12. 

0 

0 

.00 

.000 

152   24. 

0 

0 

.00 

.000 

152  48. 

0 

0 

.00 

.000 

152   72. 

0 

0 

.00 

.000 

152  120. 

0 

0 

.00 

.000 

152  192. 

0 

0 

.00 

.000 

:es    RATI 

HOURS 

VOLUME 

.0010 

.78 

.023 

.0016 

.76 

.023 

.0025 

.74 

.022 

.0040 

.71 

.021 

.0063 

.65 

.019 

.0100 

.57 

.017 

.0160 

.48 

.014 

.0250 

.39 

.010 

.0400 

.30 

.005 

.0630 

.10 

.000 

.1000 

.00 

.000 

.1600 

.00 

.000 

.2500 

.00 

.000 

.4000 

.00 

.000 

.6300 

.00 

.000 

1.0000 

.00 

.000 

1.6000 

.00 

.000 

2 .5000 

.00 

.000 

6.3000 

.00 

.000 

10.0000 

.00 

.000 

Figure  8. — Example  of  output  obtained  from  the  input  shown  in  figure  2, 
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APPENDIX  1 


FLOW  CHART  OF  THE  RUNOFF  PROGRAM 


0 


START 


(READ:LDAY(M)  ,LEAP(M)  , 
V   CLASS (J), P (J) 


_^_ 


READ:ID,NA,NB,NC, 

NQ,NP,NAF,QI,PI, 

AI,GDE,PDE,ACRE 


;'READ:Q(N)  ,\ 
\N=1,NQ    / 

^ PRINT:  \ 

I  heading; 


1 

LINE=0 

"m=° 

MF=1 

IS(*=1 

SHIFT=0 

fV0LD(J)=0 
Set\.DUR(J)=0 

4V=0 

SUM=0 

KYRT=0 

I7D=0 

QP=0 

EATE=0 

'1P> 


A> 


READ :  PC  (N),1 
N=1,NP 


(200) 

i-<dTc>^- 


--<£) 


READ:AREA(N) 
N=1,NAF 


'READ:  ID,  M,  KDAY, 

KYR,  (TH(I),  TM(I), 
,  H(I),  1=1,4) 


-® 


READ:ID,K,KDAY, 
KYR, (TH(I) ,TM(I) , 
H(I),D(I),I=1,3) 


"D0"I=1,4| 
H(I)=H(I)+SHIPT -I 


213) <   M:5 

|SHIFT=H(1)|  |h(I)=H(I)+SHIPt| 1 


Begin  Maximum  Volume  Search 
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APPENDIX  1--CONTINUED 
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APPENDIX  1--CONTINUED 


DUR(J)=DUR(J)+^t1 
VDLD  (J)  =V0LD(J)+(QOK30qiI»)*»  t  x /2- Cft  t  i*CLASS  ( J)  ) 


DnR(J)-DITR(J)+(QC0RR-CLASS(J))*At1/(()a)Bi-(30) 
VDLD(J)=VOIJ)(J)+(QCORR-CLASS(J))**2*<it1/(tQODRR-qO)*2) 


DOR(J)-Dint(J)-K<}0-ClASS(J))*atL/(QO-QOORK) 
VDI*(J)-VOLD(J)+(QO-CIJLSS(J))**2*<Jt1/(tQO-QOORK)*2) 


[CONTINUE  |— 
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APPENDIX  1-CONTINUED 


0 


"TV=L-1 


"W 


SUM+aV 


SUM=VLV 

IVp^+AV 
STORM=STORM+.AV 
QC0R2=QC0KR*1 .008*ACRE 


f      PRINT:  \ 

V^DATA  HEADING^/ 


Begin  Maximum  Volume  Search 


KDAYT=KDAY 
ZVff=0 


Ad  -cA"l 
HCk=HC2 


DFT=DF 
HCO=HCl 
QO=QCORR 


JA=3 


1:4 


-®- 


JA=1 


5r 


HRT-TH(I) 

TMIN=TM(I) 

I^I/t-1 


(921 


"© 
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APPENDIX  1--CONTINUED 


APPENDIX  2 
DIVISION  OF  THE  PROGRAM  BY  FUNCTIONS 

Division 

number6  Program  functions 

1.  Initializing  and  input. 

2.  Converting  time  and  making  preliminary  decisions. 

3.  Determining  whether  year  is  leap  year. 

4.  Computing  time  increments  and  indexing  stage  readings. 

5.  Computing  pondage  rate. 

6.  Computing  observed  discharge. 

7.  Combining  pondage  and  observed  rates. 

8.  Correcting  for  variable  drainage  area. 

9.  Computing  flow  volume. 

10.  Determining  flow  durations  and  volumes. 

11.  Accumulating  volumes. 

12.  Printing  answers. 

13.  Re-indexing  and  returning  to: 

1)  read  a  new  card, 

2)  next  value  on  last  card  read,  or 

3)  maximum  volume  search. 

14.  Maximum  volume  search. 


""Corresponding  to  divisions  on  the  flow  chart,  appendix  1, 
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FORTRAN  PROGRAM  LISTING 


PROGRAM  RUNOFF 

0DIMENSI9N    Q(200),     PC(200),    LDAY(12), 

LEAP(12),     TH(4),     TM(4),     H(4),    D(4), 

1T(1200),      SUMV(1200),     V(9),     T<?(9), 

KTD(9),  TD(9),  TAD(9),  P(12), 

2CLASS(20),       DUR(20),      V$LD(20), 

AREA(50) 

3  F<?RMAT(9F5.0) 

4  FQIRMAT     (1214) 

5  FQRMAT  (15,   312,    314,    2F4.2,    F4.1, 

2F7.2,  F7.0) 

6  FORMAT  (10F8.3) 

8  FORMAT     (15,    313,       3(2F4.0,      F6.2, 

F5.2)) 

7  FORMAT  (49H  NEGATIVE  DELT  TIME 

INCURRED  AFTER  NEXT  READ- 
ING) 
30  FORMAT  (F9.3,  F8.4) 

9  FORMAT  (15,  313,  4  (2F4.0,  F6.2)) 
153  FORMAT  (23X,  F8.4,  F10.2,  F10.3) 

15  FORMAT  (20X,     15,     F5.0,    213,    F6.2, 
F7.3) 
173  FORMAT  (27H  =ANNUAL  FL$IW  SUM- 
MARY ID=  14,  8H  YEAR  =13,  ///) 
163  FORMAT  (53H    DURATION    g»F    FL^W 

RATES  RATE  HQURS  VQLUME) 
162  FORMAT  (51H    MAXIMUM    FL$WS    ID 
INTER  M<?  DAY  H<?UR  VOLUME) 
1670  FORMAT  (113H  =  ID  M<?  DAY  YR  H0UR 
MIN      DELTTIME     HEAD     DHDT 
PQNDCG-RQ9  1B(CFS)  CFS  Q(IPH) 
STORM     DAILY     MONTHLY 
ANNUAL) 
151  FORMAT  (10F8.4) 
MA=1 
11  FORMAT  (15,     313,   2F4.0,   F9.4,   F6.2, 
F7.4,  3F8.2,  5F9.4) 
READ  4,  (LDAY(M),  M  =  1,  12) 
READ  4,  (LEAP(M),  M=  1,  12) 
READ  151,  (CLASS(J),  J=  1,  20) 
READ  3,  (P(J),  J=l,  9) 
14  READ  5,     ID,    NA,    NB,    NC,    NQ,    NP, 
NAF,  QI,  PI,  AI,  GDE,  PDE,  ACRE 
IF(ID)  160,  22,  160 
160  READ  6,  (Q(N),  N=  1,  NQ) 
PRINT  167 


LINE  =  0 
MF  =  1 
SHIFT=  0. 
DELV=  0. 
SUM=0. 
SUMVD  =  0. 
SUMVM  =  0. 

isw  =  r 

D<?  154  J=l,  20 
V$LD(J)=0. 
154DUR(J)=0. 
KYRT  =  0. 
QP=0. 
RATE  =  0. 

G0  T<?  (18,  41),  NA 
18  READ  6,  (PC(N),  N=l,  NP) 
G0  T0  (200,  41),  NB 

200  Gty  T0  (201,  20),  NC 

201  READ  6,  (AREA(N),  N=  1,  NAF) 
G$  T<?  20 

41  READ    9,    ID,    M,    KDAY  KYR,   (TH(I), 
TM(I),  H(I),  1=1,  4) 
IF  (M-55)  215,  216,  215 

215  D<?  217  1  =  1,  4 
217  H(I)=H(I)+SHIFT 

Gty  Tty  21 

216  SHIFT  =H(1) 
G<?  T<?  41 

20  READ    8,    ID,   M,   KDAY,   KYR,   (TH(I), 

TM(I),  H(I),   D(I),   1=  1,   3) 
IF(M-55)  212,  213,  212 

212  D0  214  1  =  1,  3 
214  H(I)=H(I)+SHIFT 

G<?  T<?  21 

213  SHIFT  =  H(1) 
G<?  T<?  20 

21  IF  (ID)  23,  22,  23 

22  ST0P 

23  IF  (M-99)  25,  27,  25 
27  MA  =  1 

NL  =  L-2 

MR  =  1 

MG=MF 

G#  T<?  140 
25  1=1 
48  IF(TH(I))  303,  26,  26 


23 


APPENDIX  3--CONTINUED 


303  ISW  =2 

TH(I)  =  -TH(I) 
26  HRS  =  TH(I>f(TM(I)/60.) 

IF  (HRS)  31,  28,  31 
28  1  =  1 

GO  TO  (20,  41),  NB 

31  IF(I-l)  91,  32,  91 

32  IF(KYR-KYRT)  85,  33,  85 

33  ME  =  2 

Cg>  TO  (35,  36),  MF 

35  DAY  =  KDAY+LDAY(M) 
GO  TO  37 

36  DAY  =  KDAY+LEAP(M) 

37  GO  TO  (38,  91),  MA 

38  KYRT=KYR 
MQT  =  M 
KDAYT=  KDAY 
L  =  l 

91  T(L)=(DAY*24.)+HRS 

39  GO  TO  (89,  90),  NB 

89  DF  =  D(I)/12.0 

90  HC  =  H(I) 

GO  TO  (43,  49),  ME 

43  G<?  T9  (209,  210),  MA 

209  DELT1=  0.01667 
MA  =  2 

ISW=  3 
GO  TO  47 

210  G<?  T0  (44,  45),  MG 

44  DELT2  =  T(L)+8760.-T(L-l) 
GO  TO  47 

45  DELT2  =  T(L)+8784.-T(L-l) 
47  IF(KYR-KYRT)  50,  46,  50 

46  ME  =  2 
DELD=0. 
DFT  =  DF 
HC1=  HC 
HCO=HC 
KDAYT  =  KDAY 
MOT  =  M 
Q9=0. 
TMIN=TM(I) 
HRT=TH(I) 
1=1+1 

L=L+1 
GO  TO  26 
49  DELT2  =  T(L)-T(L-1) 
IF(DELT2)  400,  400,  50 
400  PRINT  7 

GO  TO  (401,  402),  NB 
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401  IF  (I- 3)  403,  20,  20 

402  IF(I-4)  403,  41,  41 

403  I  =  1+1 
GO  TO  48 

50  HC2  =  HC 

GO  TO  (52,  57),  NA 
52  DELD2  =  DF-DFT 

IF(DELD2)  54,  220,  220 

54  DELD2=0. 

220  GO  TO  (184,  203,  204),  ISW 

203  ISW  =  3 
QOB  =  0. 
QP  =  0. 
RATE  =  0. 
GO  TO  63 

204  QP=  0. 
ISW=  1 
STORM  =  0.0 
GO  TO  94 

184  RATE  =((HCl-DELD-HCO)/DELTl+ 
(HC2-DELD2-HC1)/DELT2)/120. 
HN  =  (HC1-PDE)/PI+1.0 
N  =HN 
HCN=  N 
IF  (HCN-HN)  55,  56,  55 

55  NI  =  N+l 

QP  =  (PC(N)+(HN-HCN)*1(PC(NI)-PC(N))) 

♦RATE 
GO  TO  57 

56  QP  =  PC(N)*RATE 

57  HN  =  (HC1-GDE)/QI+1.0 
IF(HN)  94,  94,  93 

93  N  =  HN 
HCN=  N 

58  IF(HCN-HN)  59,  60,  59 

94  QOB  =  0. 
Gty  TO  61 

59  NI  =  N+l 

QOB  =Q(N>f(HN-HCN)*(Q(NI)-Q(N)) 
GO  TO  61 

60  Q^B  =  Q(N) 

61  QCORR  =  QOB+QP 
IF(QCORR)  63,  64,  64 

63  QCORR  =  0. 

64  GO  TO  (206,  205),  NC 

206  HN  =  (HC1-PDE)/AI+1.0 
N  =HN 

HCN  =N 

IF(HCN-HN)  208,  207,  208 

207  ACRE  =  AREA(N) 
GO  TO  205 
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208  NI  =  N+l 

ACRE  =  AREA(N)-KHN-HCN)*(AREA 
(NI)-AREA(N)) 
205  QC^RR  =  QC0RR/(1.OO8*ACRE) 

DELV=  (Q0+QC(/>RR)*DELTl/2. 

D9  183  J=l,  20 
175IF(QC9RR-CLASS(J))  177,  178,179 

177  IF(Q9-CLASS(J))  75,  75,  180 

180  DUR(J)  =DUR(J)+(QQ-CLASS(J))  *  DELT1/ 

(Q9-QC^RR) 
V0LD(J)  =  V0LD(J>f(Q<?-CLASS(J))* 

(Q  9-C  LAS  S(J))*D  E  L  T1/((Q  <?- 

QC0RR)*2.) 
G<?  Tty  183 

178  JB  =  1 

G<?  T<?  176 

179  JB      2 

176  IF(Q9-CLASS(J))  181,  152,  152 

181  Gty  T0  (75,  185),  JB 

185  DUR(J)  =  D  U  R(J)+(Q  C  0  R  R-CLASS(J)) 
*DELT1/(QC0RR-Q<?) 
OVQLD(J)  =  V0LD(J)+(QC0RR-CLASS 
(J))*  (QC0RR-C  LAS  S(J))*DELT1/ 
((QC0RR-Q01)*2.) 
Gq>  T0  183 
152  DUR(J)  =DUR(JHDELT1 

V0LD(J)  =  V9LD(J)+(Q#+QC<?RR)* 
DELTl/2.-DELTl*CLASS(J) 
183  CONTINUE 
75  LV=  L-l 

SUMV(LV  =  SUM  +  DELV 
SUM=SUMV(LV) 
SUMVM=  SUMVM  +  DELV 
SUMVD  =SUMVD  +  DELV 
STQRM  =  ST^RM  +  DELV 
QCQR2  =  QCQRR  *  1.008  *  ACRE 
760  PRINT  11,    ID,    M0T,    KDAYT,   KYRT, 
HRT,   TMIN,  DELT1,  HC1,  RATE, 
QP,  Q^B,  QC9R2,  1QC0RR,  ST(?RM, 
SUMVD,  SUMVM,  SUM 
LINE  =  LINE  +  1 
IF  (56-LINE)  168,  168,  169 
169  G0  Tg>  (69,  66),  ME 
168  PRINT  167 
LINE  =  0 
G0  T<?  169 
69  NL=  L-l 
MR=  2 
G<?  T0  140 


62  KYRT=KYR 
ME=2 
L=l 

T(L)=T(NL  +  2) 
SUM=0. 

67  SUMVM  =  0. 
M<?T=M 

68  KDAYT=KDAY 
SUMVD  =0. 

G0  T<?  74 
66  IF(MOT-M)  67,  71,  67 
71  IF  (KDAY-KDAYT)  68,  74,  68 
74  DELT1 =DELT2 

DFT=DF 

DELD=DELD2 

HCO=HCl 

HC1  =  HC2 

Q0  =  QC<?RR 

G<?  T0  (79,  80),  NE 

79  IF  (1-3)  81,  82,  81 

80  IF  (1-4)  81,  83,  81 

81  JA=1 
77HRT=TH(I) 

TMIN=TM(I) 

L  =  L+1 

Gq>  T<?  (17,  24,  92),  JA 

82  JA=2 
G$  Tty  11 

83  JA=3 

Gq>  Tty  11 
17  1=1  +  1 

Gq>  T0  48 
24  1=1 

G^  Tty  20 
92  1=1 

G0  T0  41 
85  ME=1 

MG  =  MF 

YR=KYR 

TEST=YR/4. 

KTEST=TEST 

TESTA =KTEST 

IF(TEST-TESTA)  88,  87,  88 

87  MF  =  2 

DAY=  KDAY+LEAP(M) 
Gg>  T0  37 

88  MF  =  1 

DAY=  KDAY  +  LDAY(M) 
G<?  T<?  37 
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140  L  =  1 
J=  1 
96V(J)  =  0. 
TO(J)  =  0. 
LL=2 
LT=  1 
KE=  1 

97  IF  ((T(LL)-T(LT))-P(J))98,99,122 

98  LL=LL  +  1 

IF  (LL-NL-1)  97,  107,  107 

99  VOL=SUMV(LL)-SUMV(LT) 

100  IF  (VQL-V(J))  106,  101,  101 

101  V(J)  =  V0L 
TO(J)=T(LT) 

106  IF(NL-LL)  107,  107,  102 

102  LL  =  LL+1 
LT=LT+1 

IF  (T(LL)-T(LL-1)-T(LT)+T(LT-1)) 
105,  99,  104 

104  KE  =  1 

GO  TO  122 

105  KE  =  2 

GO  TO  130 

107  IF(9  -  J)  108,  109,  108 

108  J  =  J+1 
GO  TO  96 

109  D<?  110  J  =  l,  9 
TD(J)=TQ(J)/24. 
KTD(J)  =  TD(J) 

110  TAD(J)  =  KTD(J) 
PRINT  173,  ID,  KYRT 
GO  TO  (118,  112),  MG 

112  PRINT  162 

DO   115  J=  1,  9 

IF  (TD(J))  545,  545,  546 

545  MONTH  =  0 
MDAY=  0 
HRS  =  0. 
GO  TO  547 

546  D0  117  MH=  1,  12 

IF  (KTD(J)-LEAP(MH))  114,  114,117 
117  CONTINUE 
MH=13 

114  MONTH  =  MH-1 
MDAY=KTD(J)-LEAP(MH-1) 

150  HRS  =  (TD(J)-TAD(J))  *  24. 

547  PRINT  15,    ID,    P(J),    MONTH,   MDAY, 

HRS,  V(J) 

115  CONTINUE 
G0  TO  171 


118  PRINT  162 
DO  120  J=l,  9 

IF  (TD(J))  542,  542,  543 

542  MONTH  =  0  * 
MDAY=  0 
HRS  =  0. 

GO  TO  544 

543  DO  H9  MH=1,  12 

IF  (KTD(J)-LDAY(MH))  121,  121,  119 

119  CONTINUE 
MH=13 

121  MONTH  =  MH-1 

MDAY=KTD(J)-LDAY(MH-1) 
HRS=(TD(  J)-TAD(J))*  24. 

544  PRINT  15,    ID,    P(J),    MONTH,  MDAY, 

HRS,  V(J) 

120  CONTINUE 
171  PRINT  163 

DO  164  J=l,  20 

PRINT  153,  CLASS(J),DUR(J),  VOLD(J) 
V0LD(J)=  0. 
DUR(J)=  0. 
164  CONTINUE 
T(1)=T(NL) 
GO  TO  (14,  211),  MR 
211  PRINT  167 
LINE  =  0 
GO  TO  62 
127  LT=LT+1 

1220  V0L=SUMV(LL-  1)  +(P(J)-T(LL-1)+ 
T(LT))*(SUMV(L  L)-SUMV(LL-l))/ 
(T(LL)-1T(LL-1))-SUMV(LT) 

123  IF  (VOL-V(J))  125,  124,  124 

124  V(J)  =  VOL 
T0(J)=T(LT) 

125  GO  TO  (126,  134),  KE 

126  REMANL=T(LL)-T(LT)-P(J) 

IF  (REMANL-T(LT+1)+T(LT))  129,  128, 
127 

128  LT=LT+1 

GO  TO  99 

129  LT=LT+1 

KE  2 
1300  V0L=SUMV(L  L)-SUMV(L  T-1)-(T(LL)- 
T(LT-1)-P(J))*(S  U  M  V(L  T)-SUMV 
(LT-1))1/(T(LT)-T(LT-1)) 

131  IF  (VOL-V(J))  125,  132,  132 

132  V(J)=V0L 
T0(J)=T(LL)-P(J) 
G9  TO  125 
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134  REMANT=P(J)-T(LL)+T(LT)  138  LL  =  LL+1 

IF  (NL-LL)  136,  107,  136  G<?  T$  99 

136  IF  (REMANT-T(LL+1)+T(LL))  139,  138,  139  LL   =  LL+1 

137  KE   =  1 

137  LL  =  LL  +  1  G<?  T<?  122 
KE   =  2  END 

G0  T<?  130  END 
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DEFINITION  OF  PROGRAM  SYMBOLS 


Flow  chart 7 
symbol 


FORTRAN 
symbol 


Definition 


AV 
Ad 
Ad„ 

Atl 
At. 


ACRE 
AI 

AREA(N) 

CLASS  (J) 

DAY 
DELV 

DELD 

DELD2 

DELT1 
DELT2 
DF 

DFT 

D(I) 

DUR(J) 
GDE 


Drainage  area  in  acres. 

Stage  interval  in  feet  between  items  in  the  variable 
drainage  area  correction  table. 

Area  of  the  watershed  in  acres  tabulated  at  uniform 
stage  increments  AI,  beginning  at  PDE. 

Runoff  rate  in  inches  per  hour.  Defines  20  points  on 
the  duration  curve.  Suggested  values  are  shown  in 
table  1. 

Day  of  year. 

Differential  volume  of  runoff  in  inches  in  time  period 
Atr 

Incremental  depth  of  precipitation  falling  in  time 
period  At  . 

Incremental  depth  of  precipitation  falling  in  time 
increment  At„. 

Time  increment  in  hours.  (See  fig.  5.) 

Time  increment  in  hours.  (See  fig.  5.) 

Accumulated  depth  of  precipitation  in  feet.  DF  = 
D(I)/12. 

The  (i-l)—  depth  of  accumulated  rainfall  in  feet. 

Accumulated  depth  of  precipitation  recorded  by  a 
nearby  raingage  (inches). 

The  time  in  hours  that  the  discharge  rate  is  equal  to 
or  greater  than  the  value  specified  by  CLASS(J). 

Gage  datum  elevation.  Elevation  of  zero  discharge 
and  first  item  in  rating  table. 


HRT 


If  different  from  FORTRAN  symbol. 


Time  in  hours  of  previous  reading. 
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Flow  chart 

FORTRAN 

symbol 

symbol 

H(I) 

HRS 

HC 

HC0 

HCO 

HC1 

HC1 

HC2 

HC2 

HCN 

HN 

Definition 


Gage  height  reading. 

Hours  of  the  day  of  the  current  gage  height  reading. 

Equivalent  to  H(I). 

Gage  height  reading.  (See  fig.  5.) 

Gage  height  reading.  (See  fig.  5.) 

Gage  height  reading.  (See  fig.  5.) 

Decimal  equivalent  of  N. 

Position  in  discharge  rating  table,  pondage  correction 
table,  or  area  correction  table. 

I  Integer  variable  indexing  TH(I),  TM(I),  H(I),  and 

D(I). 

ID  Identification  number  of  station. 

IDT  Identification  number  of  previous  card. 

ISW  Code  used  in  branching  instructions.  Normal  value 

is  1.  At  coded  beginning  of  event  (TH(I)  is  negative), 
ISW  is  set  to  2.  Its  value  then  progresses  to  3  and  is 
reset  to  1. 

J  Integer  variable  used  as  an  index. 

JA  Code  used  in  branching  instruction.  May  have  values 

of  1,  2,  or  3. 

JB  Code  used  in  branching  instruction.  May  have  values 

of  1  or  2. 

KDAY  A  two-digit  integer  designating  the  day  of  reading. 

KDAYT  Day  of  previous  card. 

KYR  Year  of  readings  (last  two  digits). 

KYRT  The  year  of  previous  card. 

KE  Integer  variable  with  possible  values  of  1  or  2  for 

branching. 
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Flow  chart  FORTRAN 

symbol  symbol  Definition 

KTEST  Integer  which  equated  to  YR  drops  off  the  decimal 

portion  of  YR. 

KTD(J)  The  day  of  year  associated  with  the  beginning  of  a 

period  of  maximum  volume,  V(J). 

L  Integer  index  for  storing  cumulative  volumes  and 

times. 

LDAY(M)  The  day  of  the  year  of  the  first  day  of  month,  M,  in  a 

non-leap  year. 

LEAP(M)  The  day  of  the  year  of  the  first  day  of  month,  M,  in  a 

leap  year. 

LINE  Count  of  lines  of  print  on  each  output  sheet. 

LL  Integer  variable  indicating  the  At^  in  which  the 

leading  edge  of  search  mesh  is  located. 

LT  Integer  variable  indicating  the  At,  in  which  the 

trailing  edge  of  search  mesh  is  located. 

LV  Index  equal  to  L-l. 

M  A  two-digit  integer  designating  the  month  of 

reading. 

MA  Integer  variable  having  possible  values  of  1  or  2.  If 

1,  card  is  first  for  a  particular  station. 

MDAY  Day  of  month  of  beginning  of  period  of  maximum 

volume,  V(J). 

ME  Integer  variable  having  possible  values  of  1  or  2.  If 

1,  reading  is  first  of  new  year. 

MF  Integer  variable  having  possible  values  of  1  or  2.  If 

2,  current  year  is  a  leap  year. 

MG  Integer  variable  having  possible  values  of  1  or  2.  If 

2,  previous  year  was  a  leap  year. 

MH  Integer  index;  refers  to  month  of  year. 

M^ZNTH  Month  of  beginning  of  period  of  maximum  volume, 

V(J). 

M0T  Month  of  previous  card. 
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Flow  chart  FORTRAN 

symbol  symbol  Definition 

MR  Integer  variable  having  possible  values  of  1  or  2. 

Used  in  branching  instruction. 

N  Integer  equivalent  of  HN  (truncates  decimal  portion 

of  HN).  Lower  limit  of  interpolation  interval  in  rating 
tables. 

NA  A  one-digit  integer  having  possible  values  of  1  or  2. 

If  1,  pondage  correction  is  to  be  applied. 

NAF  A  three-digit  integer  designating  the  number  of  items 

in  the  variable  drainage  area  correction  table. 

NB  A  one-digit  integer  having  possible  values  of  1  or  2.  If 

1,  correction  is  made  for  rainfall-on-the-pond. 

NC  A  one-digit  integer  having  possible  values  of  1  or  2.  If 

1,  correction  is  made  for  varying  drainage  area  as 
pond  stage  changes. 

NI  Upper  limit  of  interpolation  interval  in  rating  tables. 

NL  Terminating  value  of  L. 

NP  A  three-digit  integer  designating  the  number  of  items 

in  the  pondage  correction  table. 

NQ  A  three-digit  integer  designating  the  number  of  items 

in  the  discharge  rating  table. 

PC(N)  Pondage  correction  in  cfs  for  a  rate  of  change  in 

stage  of  one  foot  per  minute. 

PC(NI)  Same  as  PC(N),  except  that  NI  =  N+l. 

PDE  Pond  datum  elevation.  Elevation  of  first  item  in 

pondage  correction  table. 

PI  Stage  interval  in  feet  between  items  in  the  pondage 

correction  table. 

P(J)  The  time  period  in  hours  for  which  maximum  volumes 

are  computed.  (See  table  3.) 

^0  Q0  Discharge  rate  in  cfs  corresponding  to  HCO. 

QcjZJRR  QC0RR  Discharge  rate  corrected  for  pondage. 
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Flow  chart 
symbol 


FORTRAN 
symbol 


Definition 


V, 


D 


V. 


M 


QCOR2 

QI 

Q(N) 

Q(NI) 

Q0B 

QP 

RATE 
REMANL 

REMANT 

SHIFT 

STpRM 
SUM 
SUMVD 
SUMV(L) 

SUMVM 

TAD(J) 

TD(J) 

TEST 

TESTA 

TH(I) 


Discharge  rate  corrected  for  pondage,  cubic  feet  per 
second 

Stage  interval  in  feet  between  items  in  the  discharge 
rating  table. 

Discharge  values  in  cfs  in  rating  table. 

Same  as  Q(N)  except  that  NI=  N+l. 

Discharge  rate  in  cfs  for  a  specified  head  reading. 
Not  corrected  for  pondage. 

Pondage  correction  in  cfs  for  a  computed  rate  of 
change  in  stage. 

Rate  of  change  in  stage,  feet  per  minute 

Defined  as: 

REMANL  =  T(LL)  -  T(LT)  -  P(J). 

Defined  as; 

REMANT  =  P(J)  -  T(LL)  +  T(LT). 

A  correction  to  be  added  algebraically  to  all  head 
readings. 

Accumulated  storm  flow  in  inches. 

Accumulated  annual  flow  volume  in  inches. 

Accumulated  daily  flow  volume  in  inches. 

Accumulated  flow  in  inches,  indexed  and  stored  for 
maximum  volume  search. 

Accumulated  monthly  flow  volume  in  inches. 

Decimal  equivalent  of  KTD(J). 

Defined  as: 

TD(J)  =  TO(J)/24. 

YR  divided  by  4  (used  in  test  for  leap  year). 

Decimal  variable  corresponding  to  KTEST. 

Military  time  in  hours  of  i     reading. 
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APPENDIX  4--CONTINUED 


Flow  chart 
symbol 


FORTRAN 
symbol 


Definition 


T(L) 

T(LL) 

T(LT) 

TM(I) 
TMIN 
T(NL) 
Tp(J) 

V(J) 

VOL 
VOLD(J) 

YR 


Time,  in  hours  of  the  year,  associated  with  accumu- 
lative flow  values. 

Time,  in  hours  of  year,  associated  with  the  leading 
edge  of  the  search  mesh. 

Time,  in  hours  of  year,  associated  with  the  trailing 
edge  of  the  search  mesh. 

Time  in  minutes  of  i —  reading. 

Time  in  minutes  of  previous  reading. 

Same  as  T(L).  NL  is  the  terminating  value  of  L. 

Time  in  hours  of  year  of  beginning  of  period  of 
maximum  volume,  V(J). 

Maximum  flow  volume  occurring  in  time  interval 
P(J). 

A  flow  volume  occurring  in  time  interval  P(J). 

Volume  of  flow  in  inches  occurring  at  rates  = 
CLASS  (J). 

Decimal  variable  corresponding  to  KYR. 
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